Aging in a Chaotic System 
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Abstract 

We demonstrate aging behavior in a simple non-linear sys- 
tem. Our model is a chaotic map which generates determin- 
istically sub-diffusion. Asymptotic behaviors of the diffusion 
process are described using aging continuous time random 
walks, introduced previously to model diffusion in glasses. 
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Aging behavior is found in complex dynamical systems 
like spin glasses, glasses, and polymers These systems 
include many interacting sub-units and are disordered. 
In this letter we will demonstrate that aging can be found 
also in low dimensional non-linear systems. Specifically 
we will show that deterministic diffusion generated by 
one dimensional maps exhibits aging. And that statis- 
tical properties of the corresponding trajectories can be 
analyzed using aging continuous time random walks (AC- 
TRW), introduced by Monthus and Bouchaud Q in the 
context of aging dynamics in glasses (see also (j| and Ref. 
therein) . 

There exist several methods to investigate aging. One 
method is to start a dynamical process at time t = — t a , 
then at time t — add a small perturbation to the sys- 
tem. One eventually measures the response at some time 
t > 0. Alternatively one can measure correlation func- 
tions between physical quantities at time t = and time 
t, after aging the system in the interval (— t a , 0). We use 
the latter approach. A system exhibits aging if its dy- 
namical properties depend on t and t a even in the limit 
when both are long. Of-course many systems do not 
exhibit aging, namely when t > r, where t is a charac- 
teristic time scale of the problem, dynamical properties 
of the system are independent of the aging time ta- 
in, many cases trajectories generated by determinis- 
tic systems, such as low dimensional Hamiltonians or 
maps, are highly irregular and noisy. For an observer 
these trajectories seem to be generated by a stochastic; 
rather than a deterministic mechanism. Hence analysis 
of chaotic trajectories generated deterministically is often 
based on random walk concepts jl] |i) . It is well known |]] 
that both conservative and dissipative deterministic sys- 
tems may generate normal (a = 1) or anomalous (a^l) 
diffusion for a coordinate x 

(x 2 )~e, (i) 

where the average is over a set of initial conditions (see 



details below). The anomalous behavior is due for exam- 
ple to long trapping events in vicinity of unstable fixed 
points or stickiness near islands in phase space [|l0| . 
Here we shall investigate a deterministic system which 
generates sub diffusion behavior a < 1 and show how 
rich dynamical features emerge when aging is included. 
More generally, we demonstrate that aging can be used 
to probe dynamics in deterministic systems. 

Probably the simplest theoretical tool which gener- 
ates normal and anomalous diffusion are one dimensional 
maps 



x t+ i =x t + F(x t ) 



(2) 



with the following symmetry properties of F(x): (i) F(x) 
is periodic with a periodicity interval set to 1, F(x) = 
F(x + N), where N is an integer, (m) F(x) has inversion 
anti-symmetry; namely, F (x) — —F(—x). The study of 
these maps was motivated by the assumption that they 
capture essential features of a driven damped motion in a 
periodic potential [|ll) . Geisel and Thomae ||] considered 
a rather wide family of such maps which behave as 



F(x) 



for x 



-0, 



(3) 



where z > 1. Variations of these maps have been in- 
vestigated by taking into account: time dependent noise 
Jl2[ , quenched disorder |l3| , and additional uniform bias 
which breaks the symmetry of the map Jl4| . We use the 
map 



F(x) = (2xY 



< x < 



(4) 



which together with the symmetry properties of the map 
define the mapping for all x. In Fig. [l] we show the map 
for three unit cells. 

To investigate aging, e.g. numerically, we choose an 
initial condition x_ ta which is chosen randomly and uni- 
formly in the interval —1/2 < x^ ta < 1/2. The quan- 
tity of interest is the displacement in the interval (0, t), 
x = Xt — xq which is obtained using the map Eq. (|^). 
Previous work ]§|,|T5t considered the non-aging regime, 
namely t a = 0. 

In an ongoing process a walker following the itera- 
tion rules may get stuck close to the vicinity of unsta- 
ble fixed points of the map (see Fig. |l|). It has been 
shown, both analytically and numerically, that probabil- 
ity density function (PDF) of escape times of trajecto- 
ries from the vicinity of the fixed points decays like a 



1 



power law |g. To see this one considers the dynamics 
in half a unit cell, say < x < 1/2. Assume that at 
time t = the particle is on x* residing in vicinity of 
the the fixed point x — 0. Close to the fixed point we 
may approximate the map Eq. (|2|) with the differential 
equation dx/dt — F(x). Hence the escape time from x* 
to a boundary on b (x* < b < 1/2) is 



t 



b dx 



using Eq. (§) 



1 

t ~ - 

a 



-z+l 



b -z + l 



1 



1 



(•5) 



(6) 



The PDF of escape times tp(t) is related to the unknown 
PDF of injection points T](x*), through the chain rule 
if}(t) = r)(x*)\dx* /dt\. Expanding r)(x*) around the un- 
stable fixed point x* = one finds that for large escape 
times 



T(-a) 



-r 



(z-iy 



(7) 



where A depends on the PDF of injection points. Note 
that when z > 2 corresponding to a < 1, the average 
escape time diverges. 




later using numerical simulations. As mentioned start of 
the ACTRW process is at t = —t a and our goal is to find 
the ACTRW Green function P(x, t a , t). 

In ACTRW we must introduce the distribution of the 
first waiting time t\: the time elapsing between start 
of observation at t = and the first jump event in the 
interval (0,t). Let h ta (h) be the PDF of t x . Let h s (u) 
be the double Laplace transform of /it a (ii) 



h s (u) = / dh / dt a ht a (h) e 
Jo Jo 



(8) 



then according to theory of fractal renewal processes |T(| 
1 ip(s) — ip(u) 



h s (u) = 

1 - ip(s) u — s 

When a < 1 in Eq. (0) (i.e., z > 2) 

sin(7ra) t" 



ht a (h) 



TT t? (*l + t a ) ' 



(9) 



(10) 



which is valid in the long aging time limit t a >> A 1 /". 
Note that Eq. ([l(]) is independent of the exact form of 
ip(t), besides the exponent a. When a — * 1 the mass of 
the PDF ht a (t±) is concentrated in the vicinity of t\ — » 0, 
as expected from a 'normal process'. In what follows we 
will also use the double Laplace transform of Eq. (|1C 



h s (u) 



s a (u — s) 



(11) 



We have checked numerically the predictions of Eq. 
( |l0| ) for z = 3, analyzing trajectories generated by the 
map Eq. (^) with three different aging times. In Fig. ^ 
we show the probability of making at-least one step in 
the interval (Q,t): J Q ht a (t)dt = 1 —po(t a ,t). The results 
show a good agreement between numerical results and 
the theoretical prediction Eq. ([To]) without fitting. Fig. 
§ clearly demonstrates that as the aging time becomes 
larger the time for the first jumping event, from one cell 
to its neighbor, becomes larger in statistical sense. 



FIG. 1. The map x t+ i = x% + F(x t ), defined by Eq. (g) 
with z — 3. The linear dash-dot curve is xt+i = Xt. The 
unstable fixed points are on xt — 0, 1, 2. 

To consider stochastic properties of the aging dynamics 
we investigate now ACTRW, deriving an explicit expres- 
sion for the asymptotic behavior of the Green function. 
Specifically we consider a one-dimensional nearest neigh- 
bor random walk. Each lattice point corresponds to a 
cell of the iterated maps. Waiting times on each lattice 
point are assumed to be described by ip(t). Note that 
after each jumping event we assume that the process is 
renewed, namely, we neglect correlation between motions 
in neighboring cells. This assumption will be justified 
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FIG. 2. We show the probability of making at least one 
step in a time interval (0, t) for different aging times specified 
in the figure. The solid curve is the theoretical prediction Eq. 
(fi"o|), the dotted, dashed, and dot dashed curves are obtained 
from numerical solution of the map with 2 = 3. 

We now investigate the ACTRW Green function. Let 
Pn{t a , t) be the probability of making n steps in the time 
interval (0,i). Let P(k,s,u) be the double-Laplace - 
Fourier transform (x — > k, t a — > s, t — > u) of P(x,t a ,t), 
then 



Pn{s,u) COS" (k) 



(12) 




where p n (s,u) is the double Laplace transform of 
Pn(ta,t)- In Eq. ( jO ) cos"(fc) is the characteristic func- 
tion of a random walk with exactly n steps. Using con- 
volution theorem of Laplace transform 



sh s (u) 



Pn(s,U) 



i,(u) 



Q 



n > 1. 



(13) 



Hence inserting Eq. (|T^) in Eq. (|l2|), using Eq. (g) and 
summing we find 



P(k, s, u) 



1 | [^(M)-V»(a)][l-cos(fc)] 
su u (u — s) [1 — ^ (s)] [1 — "0 ( M ) cos (&)] 

(14) 



We note that only if the underlying process is a Poisson 
process, the Green function P(x,t a ,t) is independent of 

t a - 

By differentiating Eq. (|l4|) with respect to k twice and 
setting k — 0, we obtain the mean square displacement 
of the random walk 



(x 2 (s,u)) 



h s (u) 



u [1 — 1p (it)] 



From Eq. (j^) we have ip(u) = 1 — Au a ■ • •, when u is 
small. Using Tauberian theorem one can show that for 
t, t a > A 1 /" 



(X 2 (t«,t)> 



1 



1 



AT(l + a) 



[(* + *-)"-*?] 



(16) 



For times t >> t a we recover standard CTRW behavior, 
{x 2 (t a ,t)) oc t a @. For t/i a << 1 we find (x 2 (t,t a )) cx 
t/t a ~ a , hence as t a becomes larger the diffusion in this 
regime is slowed down. 

This ACTRW behavior is shown in Fig. || for the iter- 
ated map. We show (x 2 (t a ,t)) versus time for different 
values of t a . Good agreement between ACTRW and the 
numerical simulations is found. All the curves in Fig. 
H converge into (x 2 (t a ,t)) ~ t a for long forward times 
ij while at shorter times the diffusion is clearly slowed 
down as the aging time is increased. 



FIG. 3. We show the mean square displacement versus the 
forward time t for different aging times t a , ta — triangle, 
t a = 1000 diamond, t a = 10000 circle, and t a = 60000 star. 
The solid curves are the theoretical prediction Eq. ( |l6| ) . Here 

2 = 3. 

To investigate properties of the Green function 
P(x,t a ,t) in the limit of long t and t a we consider 
the continuum approximation of Eq. ( |I4| ) (soon to be 
briefly justified). We insert the large wave length ex- 
pansion cos(fc) = 1 — fc 2 /2 and low frequency expansion 
ip(v) = l- Au a in Eq. Q 



P (k, s, u) 



') Au c 



u(u — s) 



(u - s) Au a + k 2 /2 

(17) 



Inverting to the double time (t, t a ) - real space x domain 
we find that the Green function is a sum of two terms: 



P{x,t a ,t)~ Po (t a ,t)S(x) + 



(15) sin (na) 




where in this limit 



Po (*»>*) 



sin (jra) 



dx 



t/t a x a (l + x) 



(19) 



The first term on the right hand side of Eq. ( |18| ) is a 
singular term. It corresponds to a random walk which 
did not make a jump in the time interval (0, t). The 
symbol <g> in the second term in Eq. ( [Lsl) is the Laplace 
convolution operator with respect to the forward time t, 
I a /2(f) is the one sided Levy stable PDF, whose Laplace 
pair is exp(— u a l 2 ). It is easy to see that the asymptotic 
solution Eq. ( |l8| ) is non negative and normalized [proof: 
set k = in Eq. (jTfl)]. In the limit a — > 1 we get 
a Gaussian Green function which is independent of t a , 
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hence in the normal transport regime no aging behavior 
is found [proof: set a = 1 in Eq. ( |i~7| ) 
To justify the approximation Eq. 



18|) we have calcu- 



lated the even moments of the ACTRW exactly, in u, s 
space, using Eq. ( |l4| ) [i.e., {x 2 (s,u)) Eq. (|l5|), {x 4 (s,u)), 
etc]. Then using Tauberian theorems we can show that 
for long t and long t a , the ratio t/t a being arbitrary, 
these moments and the moments obtained directly from 
the approximation Eq. ( |i~8| ) are identical. In this way we 
have justified the continuum approximation. The details 
of this calculation and generalizations to other classes of 
random walks will be published elsewhere. 

The behavior of the Green function Eq. (|i~8| ) is shown 
in Fig. (||) . Not shown is the behavior on the origin which 
exhibits a singular behavior [i.e, the 5(x) term in Eq. 
(If)]. The behavior of this singular term was displayed 
already in Fig. ^. A good agreement between simulations 
and the ACTRW Green function is obtained. 

Using Eq. (|Ig| ) we find the non-singular part of the 
ACTRW on the origin 



P(x,t a ,t)\ x=0 = t- a/2 g 



to 



(20) 



where 



Q /2 sin(Tra) t X J ..{x-y) 



-a/2 



g(x)=x a/ \ , — '— d<A a -\ . (21) 

JK ' 2irT (1 - a/2) J Q y (l + y)y a V 1 

Hence we find the asymptotic behaviors 

( t \ 1 " 
3a72j \t^J t ^t a 



P(x,t a ,t)\ x=0 



2r(a)r(2-3o 



2r(l-a/2) 



t » t a . 

(22) 



In the limit of t ^> t a we recover the standard CTRW 
behavior |l7|| . Note however that the convergence of 
ACTRW to standard CTRW behavior may become ex- 
tremely slow when a is small. For example when a = 
1/12, and t/t a — 10 s , large deviations from the asymp- 
totic limit are found: 2r(l - a/2)t a / 2 P(x,t a ,t)\ x=0 = 
0.788 instead of the asymptotic value 1. Since the limit 
a — > is important for several systems JL8| |, (i.e., log- 
arithmic diffusion as found in Sinai's model), it is clear 
that aging can be important also when t > t a . 




FIG. 4. We show the Green function obtained from 
simulation of the deterministic map in scaling form for: (i) 
t a = 10 4 ,t = 10 3 , circles, (ii) t a = W 4 ,t a = 10 4 , triangles, 
(iii) and t — 4 * 10 J ,t a = 0, diamonds. The curves are our 
theoretical results which are in good agreement with the sim- 
ulations with z — 3. 

To conclude, we have demonstrated for the first time 
aging behavior in diffusion generated by iterated maps. 
In the regime t < t a we observe a slowing down of the 
dynamics as the aging time is increased. In the limit 
t >> t a we obtain previous results, the convergence to- 
wards this regime can become extremely slow when a is 
small. The aging dynamics of the iterated maps exhibits 
a universal behavior in the sense that the Green func- 
tion, P(x, t a , t) does not depend on precise shape of 
besides the exponent a — [z — 1) _1 < 1. We believe 
that aging dynamics can be found also in other low di- 
mensional deterministic systems, an issue left for future 
research. 
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